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^ ; Abstract. 

O , In this proceeding, we present the resuhs of a semi-analytic study of CDM substructure as a 

function of the primordial power spectrum. We apply our method to several "tilted" models in the 
ACDM framework, with n ~ 1 . 1 — 0.8, or CTj^ ~ 1 .2 — 0.65 when normalized to COBE. We also study 
a more extreme, warm dark matter-like spectrum that is sharply truncated below a scale of ^ 0.3 
h^' Mpc (~ 10^" h^' Mq). Contrary to some expectations, we show that the mass fraction of halo 
substructure is not a strong function of spectral slope, so it likely will be difficult to constrain tilt 
using flux ratios of gravitationally lensed quasars. On the positive side, all of our CDM-type models 
^> ' yield projected mass fractions that are in good agreement with strong lensing estimates: / « 1 .5% at 

, M ^ 10**Mp . The truncated model produces a significantly smaller fraction, / « 0.3%, suggesting 

that warm dark matter-like spectra are disfavored and potentially may be distinguished from CDM 
spectra using lensing. We also discuss the issue of dwarf satellite abundances, with emphasis on the 
cosmological dependence of the map between the observed central velocity dispersions of Milky 
' Way satellites and the maximum circular velocities of their host halos. In agreement with earlier 

(«~^ \ work, we find that standard ACDM over-predicts the estimated count of Milky Way satellites at fixed 

V'max by an order of magnitude, but tilted models do better because subhalos are less concentrated. 
Interestingly, under the assumption that dwarfs have isotropic velocity dispersion tensors, models 
with significantly tilted primordial power spectra {e.g., n^S 0.85, (Tj^^S 0.7) may underpredict the 
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O \ number of large Milky Way satelUtes with Vmax^ 40 km s ' . 



1. INTRODUCTION 

The ACDM model of a flat Universe dominated by cold, coUisionless dark matter 
(CDM), and a cosmological constant (A) has emerged as the standard framework for 
the growth of cosmic structure. With 0.3, h 0.7, and a nearly scale-invariant 

primordial spectrum of adiabatic density perturbations {P{k) k", n ^ 1), ACDM is 
remarkably successful at reproducing large scale observations. In contrast, this paradigm 
faces several challenges on galactic and sub-galactic scales [1, 2]. In Zentner & Bullock 
[3] (ZB), we emphasized that inflation does not predict exactly scale-invariant {i.e., 
n = 1) primordial spectra. Many models of inflation predict "tilted" spectra {n ^ 1), 
spectral index "running" (dn/dlnfc ^ 0), or other deviations from scale-invariance that 
have dramatic consequences on small scales. We showed that spectra with tilts of n ~ 0.9 
and/or running and fixed by COBE on large scales can greatly reduce the predicted 
central densities of dark matter halos, alleviating the "central density problem" plaguing 
ACDM. Further, the neighborhood of Cg ~ 0.75 implied by these tilts is provocatively 
close to many recent estimates of "low" values [4]. 

In this proceeding, we report on results from follow-up work to ZB. We study the 



TABLE 1. Initial power spectra from the inflationary models discussed in ZB. 
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Scale-invariant 
Inverted Power Law 



n= 1 
IPL4 
RMI 
RMII 
RMIII 
BSI 
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~0.94 
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- -0.004 

- -0.001 
~ -0.001 
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-0.95 
~0.83 
-0.65 
-0.75 
- 1.21 
-0.97 



Running-mass model I 
Running-mass model II 
Running-mass model III 
Broken scale-invariant 
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dependence of CDM halo substructure on the primordial power spectrum (PPS). Our 
models of the PPS are the same as those in ZB. We COBE normalize all spectra 
and we assume a cosmological model with = 1 — = 0.3, Cl^h^ = 0.02, and 
h = 0.72. The important characteristics of each input spectrum are summarized in Table 
1. Numerical simulations cannot have both the resolution and the statistics needed to 
study substructure so we model substructure semi-analytically using host halo merger 
histories [5] and a scheme for approximating subhalo orbits and tidal mass loss. Our 
model expands on previous work by Bullock et al. [6] and Taylor and Babul [7]. We 
calibrated our model against available data from N-body simulations; nevertheless, our 
results must be regarded as preliminary estimates to be verified by extensive N-body 
work. We present results based on 100 merger tree realizations. We give a detailed 
description of our model and further results in a forthcoming paper [8]. 



Efforts have been made to use flux ratios in multiply-imaged quasars to detect sub- 
structure in galactic halos and to use these measurements to constrain cosmology. In 
particular, Dalai and Kochanek [9] (DK) considered bounds on the PPS. As such, it is 
important to understand the theoretical predictions for halo substructure as a function of 
the PPS and, more generally, substructure distributions and characteristics as a function 
of cosmology. 

Our results on the substructure mass fraction and the PPS are summarized in Figure 
1. DK took a typical lens mass of 3 x 10^^ Mq and the lenses in their sample have 
a median redshift of ^ 0.6, so we present results for a 3 x 10^2 M^,. halo at z = 0.6; 
however, our results do not change appreciably as a function of mass or redshift. Lensing 
measurements are sensitive to the mass fraction in substructure projected onto the 
plane of the lens at a halo-centric distance of order the Einstein radius, 7?g ~ 5 kpc. 
Consequently, we show in Fig. 1 the mass fraction in substructure for the entire halo and 
the mass fraction in substructure in a 2D projection of radius R — 10 kpc. 

Notice that the substructure mass fraction is not a strong function of tilt and/or 
running. In tilted models, host halos are less concentrated and accrete their substructure 
later, and this compensates for the fact that the substructures are more fragile, and 
more easily destroyed by tides. It will be difficult to use substructure measurements 
to constrain these parameters. Only the BSI model, with a sharp drop in power at ~ 10^° 



2. SUBSTRUCTURE MASS FRACTIONS 
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FIGURE 1. The fraction of the host halo mass bound up in substructures of mass between IO'^Mq and 
M as a function of M. (a) The mass fraction in substructure for the entire halo, (b) The mass fraction in a 
2D halo-centric cylindrical projection of radius R=IQ kpc. The lines represent the average mass fractions 
and the errorbars show the dispersion among the 100 realizations. The models are labeled in each panel. 

Mq, shows deviation from the n = \ model that is significant compared to the scatter. It 
may be possible to constrain models with such an abrupt break (e.g., warm dark matter). 
DK found the halo mass fraction bound up in substructures of mass M^^^< 10^ Mq to be 
0.006< /sat~ 0.07. All of our models are consistent with this bound, but the truncated 
model is just at the edge of the allowed region. 

3. THE DWARF SATELLITES 

The "dwarf satellite problem", namely that ACDM predicts roughly an order of mag- 
nitude more halos with Vmax^ 40 km s~^ than observed Milky Way (MW) satellites, is 
an often-discussed challenge to ACDM [1]. Stoehr et al. [10] (S02) and Hayashi et al. 
[11] (H02) proposed that substructure halos may be significantly less concentrated than 
comparable field halos due to tidal effects. This implies that the values of Vmax that cor- 
respond to the observed central velocity dispersions, o^, of the MW satellites are larger 
than the values inferred by other authors. One must be cautious. Mass redistribution in 
subhalos is quite sensitive to the mass resolution of the simulation (S02) and the veloc- 
ity function (VF) of satellite halos is sensitive to the initial concentrations and accretion 
times of substructure (H02). Our semi-analytic model represents one extreme; we do not 
allow for redistribution of mass within a subhalo's tidal limit. Using our model, we can 
also quantify the cosmology dependence of the mapping between and Vmax- We have 
assumed that CDM halos can be well described by NFW [12] profiles with a particular 
^max and Rra&x (the radius at which V^ax is attained) and calculated all combinations of 
^max and i?max that lead to the observed values of for each of the MW satellites. We 
have assumed that the stars have isotropic dispersion tensors and that the stellar distri- 
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FIGURE 2. (a) The lower group of points represent a scatter plot of Vmax vs. Rmax for 10 realizations of 
the n = 1 model. The upper points correspond to RM I. The lines show the region that yields an observed 
value of (7* = 6.8 ± 1.6km s~^ for Carina. The thick solid line corresponds to the central value of a* while 
the thin soUd Unes correspond to the 1 <y errors. Consistency demands that Carina resides in a halo that has 
structural properties that lie in the region of overlap between the thin solid lines and the scattered points 
for each cosmology, (b) The predicted VFs (lines) and scatter for the n = 1 and RM I models along with 
the "observed" VF (shapes) for each model inferred from the observed values of (7*. Squares represent 
satellite velocities that would be implied if their halo profiles reflect halos in the n = 1 model, triangles 
correspond to the RM I expectations. 



butions are given by King profiles [13] with parameters given by Mateo [14]. 

Results for Carina are shown in Figure 2, along with a scatter plot of Vmax vs. Rmax 
for the surviving satellites in 10 realizations of a M^^^^ = 1.4 x 10^^ Mq host halo at 
z = for the n = 1 and RM I models. This plot shows how allowing for less concentrated 
halos helps to alleviate the dwarf satellite problem. Less concentrated halos (larger Rmax) 
require a larger Vmax because the outer stellar radius of Carina is much smaller than the 
radius at which the halo's rotation curve peaks. Larger halos are intrinsically scarcer 
objects, helping to explain the paucity of dwarf satellites. Feedback mechanisms (e.g., 
[6]) then explain the dearth of smaller halos. Figure 2 also demonstrates that the mapping 
between Oy, and y^ax is dependent upon cosmology, in particular the PPS, and so the 
same observational data imply a cosmology dependent "observed" VF. Our estimates 
for dwarf velocities in the n = 1 case compare well to the estimates made by Klypin and 
collaborators [1]. 

In the right panel of Fig. 2 we present the predicted VFs along with separate "ob- 
served" velocity functions for the n — I and RM I models. The RM I VF is a factor 
of ~ 2 lower than the n = 1 VF mainly because typical halos are less concentrated in 
this model, so that Vmax is lower at a given mass. Also notice that the "observed" VF 
is shifted significantly higher at high Vmax- This suggests that significantly tilted power 
spectra 0.85 may underpredict the number of MW satellites at high Vmax- More- 
over, subhalos are likely less centrally concentrated than field halos (as suggested by 
S02 and H02), and this serves to make the underprediction more pronounced. However, 



we recommend circumspection. Our results concerning VFs are sensitive to several as- 
sumptions such as the isotropy of the dispersion tensor. 



REFERENCES 

1. Moore, B. etal, ApJ, 524, L19-L22 (1999); Klypin, A. A. etal, ApJ, 522, 82 (1999). 

2. Debattista, V. P. and Sellwood, J. A., ApJ, 543, 704 (2000); Keeton, C. R., ApJ, 561, 46 (2001); Alam, 
S. M. K. et al, ApJ, 572, 34 (2002); McGaugh, S. et ai, ApJ, 584, in press, (2003). 

3. Zentner, A. R. and Bullock, J. S., PRD, 66, 043003 (2002), (ZB). 

4. Viana, R T. et ai, ApJ, 569, 75; Bahcall, N. et ai, ApJ, 585, in press (2002); Schuecker, R et ai, A&A, 
in press (2002); Allen, S. W. et ai, astro-ph/0208394, (2002); Brown, M. L. et ai, astro-ph/0210213, 
(2002); Hamana, I. et ai, astro-ph/02 10450 (2002); Peirpaoli, W. J. etal., astro-ph/0210567, (2002). 

5. Somerville, R. S. and Kolatt, T. S., MNRAS, 305, 1, (1999). 

6. Bullock, J. S., Kravtsov, A. V., and Weinberg, D. H., ApJ, 539, 517, (2000). 

7. Taylor, J. E., and Babul, A., ApJ, 559, 716, (2001). 

8. Zentner, A. R. and Bullock, J. S., in preparation. 

9. Dalai, N. and Kochanek, C, ApJ, 572, 25, (2001); Dalai, N. and Kochanek, C, PRL, submitted, (2002). 

10. Stoehr, F. et ai, MNRAS, 335, L84, (2002). 

11. Hayashi, E. etal, ApJ, in press, astro-ph/0203004, (2002). 

12. Navarro, J. R, Frenk, C. S., and White, S. D. M., ApJ, 490, 493, (1997). 

13. King,!., ApJ, 67,471,(1962). 

14. Mateo, M., Ann. Rev. Astron. Astrophys., 36, 435, (1998). 



